Quantum Spin Effect and Short-Range Order above the Curie Temperature 
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Abstract 

Using quantum Heisenberg model calculations with Green's function technique generalized for arbitrary spins, we found that 
for a system of small spins the quantum spin effects significantly contribute to the magnetic short-range order and strongly 
affect physical properties of magnets. The spin dynamics investigation confirms that these quantum spin effects favor the 
persistence of propagating spin-wave excitations above the Curie temperature. Our investigation suggests a reconsideration of 
prevailing point of view on finite temperature magnetism to include quantum effects and the magnetic short-range order. 



It is a long-standing debate over the nature of the para- 
magnetic (PM) state of the ferromagnetic (FM) materi- 
als, particularly in the transition metals. Early inelastic 
neutron experiments 0, 0, 0] determined the persistence 
of spin-wave (SW) like modes above the Curie tempera- 
ture (Tc) in both Ni and Fe, and these modes were in- 
terpreted as the evidence of considerable magnetic short- 
range order (MSRO) in the PM state 3 . The presence 
of MSRO was later supported by the spin and angle- 
resolved photoemission studies 0, |(| . This is rather un- 
usual because majority of magnetism theories are based 
on the absence of SW excitations above Tc- Moreover, 
by applying the spherical model (SM)approximation to 
the Heisenberg model (HM) , Shastry Jhfl et al. concluded 
that in contrast with the experiment [2| , this model, with 
fairly long-ranged interactions, has little MSRO and no 
SW peaks above Tc in Fe. A Similar conclusion was 
reached by Monte Carlo (MC) spin dynamics simulation 
of the classical HM for the same systemQ. Very sophis- 
ticated technique^ also could not detect any traces of 
strong MSRO in Fe or Ni. Recent theoretical spin dy- 
namics studies 10], however, have demonstrated that a 
strong MSRO is a 'must have' property of the itinerant 
magnets while such excitations like SW exist above Tq 
in both localized and itinerant magnets. 

It is puzzling that the applications of HM failed to 
predict the expected MSRO because BCC Fe has rather 
good local moments 0, an d to a large extent HM 
should be valid ^| . What usually is omitted in the clas- 
sical MC simulation is the quantum nature of spin. In 
quantum SM Tc, the spin correlations and the suscep- 
tibility are proportional to S(S +1), while the classi- 
cal coefficient scales as S 2 ; so the quantum effect con- 
tributes a factor Qs = 1 + S^ 1 . For small S this Qs 
can be unreasonably large. For example, for S = 1/2 
Qs = 3 which leads to the unphysical correlation be- 
tween nearest-neighbor (NN) spins (Si ■ Sj)/S 2 > 1 in 
the case of NN coupling of the simple cubic (SC) struc- 
ture. The same problem appeared in Ref.0 where in the 
case of strong MSRO and S = 1 (Si ■ Sj)/S 2 = 1.64. To 
avoid such difficulties a pragmatical approach is to scale 
the relevant quantities by S(S +1), as was done in Ref. 
0. Then, in SM the scaled quantities are independent 
of S, so the quantitative results for MSRO and region of 
SW existence for S = 1 and 5 = 00 are the sameQ. Be- 



low we will demonstrate that while the classical HM can 
describe some degree of MSRO above Tc, QSE properly 
included strongly increases MSRO and affects its influ- 
ence on other physical properties. 



For the HM hamiltonian H 



the PM state, we use the second-order Green's func- 
tion (GF) technique[Il EE EE To calculate GF 
Gfj = {{S^;Sj)) u one applies twice the equation of 
motion and then decouples the high-order GF of forms 
((S' p S^;S+)) u and (((S+Sf-S-S+)Sr;S+)) u . For 
S = 1/2 in one-dimensional system Kondo and Yamaji 
(KY) decoupled them by using a correction parameter 
a |14|. Here we extend their method for arbitrary S by 
introducing the following decoupling scheme (for i =/= p 
and p 7^ I), 

(({SpS l — S p Sf, S'i }; Sj)) u — ► (1 — v s 8ix)Ci P Gfj — CuG 
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2(SfSf) are the spin correlations, {A, B} = (AB+BA)/2 
is the symmetric product of operators, and v s is a in- 
dependent constant which will be determined later. For 
5 = 1/2 spin operator identities require v s = 1, and Eq. 
(JTJ is reduced to the KY decoupling. 

Decoupling the high order GF in the equation of mo- 
tion with Eq. Q one can obtain the following expression 
for the dynamic susceptibility 



2E n 4JnC n (l-7„') 



(2) 



where n is the shell index, J n and C n are J$j and C^j, 
correspondingly. 7^ = z~ x ^ 5 (1 — e lq Sn ) with z n being 
the total number of sites on n— th shell and S n being sites 
on that shell. 

The SW excitation spectrum is 



1/2 



(3) 
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where D n = A" 1 J2 k J k ln Ck with jk and ck bein S the 
Fourier transforms of Jij and CV,-, correspondingly. At 
this stage v s = (2 — S)/3S is obtained by comparing 
Eq.© with the well-known result ui q — (J° — J q )S in 
the FM spin correlation limit C n = 2S 2 /3. 

From Eq.J2J and the spectral theorem, the spin cor- 
relation can be written as 

C*=£*^(l-72)— coth£« (4) 

— LO a LI 

n H 

With the requirements C n = l/NJ2 q C9 1% and C = 
1 / N J2 q cq = 25(5+ 1)/3, Eqs.Q and gj can be solved 
self-consistently. Tc is determined by = (\ — 
X H (0, 0)/2 ). To check the validity of our method, in 
Fig. 1 we compare our calculated Tc /X^and E c /E 
for the BCC structure with the accurate results obtained 
by the high-temperature-expansion (HTE) methods [Tsl 
Il9| and the SM results in the NN coupling case. Here 
T<¥ F is the Curie temperature in the mean field (MF) 
approximation, E c and E are the total energies at Tc 
and zero temperature, correspondingly. The parameter 
E c /E is a proper measure of MSRO at Tc, and in the 
NN coupling case E c /E is identical to the average cosine 
of angles between NN spins. In the MF approximation 
there is no MSRO (E c /E = 0) at and above T^ F . The 
existence of MSRO suppresses Tc with respect to T(f F . 
Such suppression exists also in the SM and is identical 
for all 5 in that case. In more accurate calculations, 
however, Tc is more suppressed at smaller S. 

The MSRO parameter E c /E$ demonstrates the in- 
crease in MSRO for smaller S. Although in the SM in- 
creases even faster (E c /Eq cx Qs), this quantity already 
is not well defined owing to the appearance of E c /En > 1 
e.g. for the SC structure for 5=1/2 and in Ref.[7| for 
5 = 1. At this stage the scaling should be introduced 
which leads to the elimination of real QSE. 

Our formalism allows to obtain the following important 
result for Tc for 5 = oo 

T c = aT c M = 3T c iF /(2F + 1), (5) 

where a = 3F/(2F + 1) with F = A^EgC 1 - 
7^)~ 1 . This new and transparent expression provides 
another immediate and accurate check of applicability 
of our generalized GF formalism. For instance, it gives 
T c /(JiS 2 ) = 1.49,2.11, and 3.25 for SC, BCC and FCC 
structures which are very close to the corresponding HTE 
results 1.45,2.06, and 3.18pj). Eq.© clearly indicates 
the importance of the correction parameter a introduced 
above in the GF decoupling. For 5* = oo the parameter 
E c /Eq = 1 — F^ 1 is the same as the one obtained in the 
SM. 

The good agreement between our and HTE results in- 
dicates the applicability of this formalism for the case 
of arbitrary S and NN interaction. We also studied a 
Heisenberg hamiltonian corresponding to a realistic ma- 
terial: we used extended (four NN) interactions in BCC 



Fe: J 2 /Jr = 0.5221, J z /J x = 0.0056, and J4/J1 = 
-0.0879(23, where J X S 2 = 2.44 mRy For S = oo 
MC simulation gives T^ C /T^ F = 0.68 - 0.70,and 
Ec/E = 0.38 - 0.41. In SM T^ M /T^ F = 0.59 and 
Ec/Eq = 0.40Q5 for all S. In our formalism, for 5 = 
1/2, 1 and 00, T C /T^ F = 0.51,0.57 and 0.67, and their 
Ec/Eq = 0.79, 0.68 and 0.41, correspondingly. At 5 = 1, 
T^ F = 2334K is more than twice higher than the exper- 
imental 1040K of Fe, our calculated Tc is suppressed to 
the much lower value 1330K. Comparing with Fig.l, one 
can see the additional suppression of Tc with Ec / Eq be- 
ing considerably larger, thus indicating stronger MSRO 
than in the corresponding NN coupling case. However, 
the parameters in Ref . [2(| have been obtained in the long- 
wavelength approximation and can only describe a small 
MSRO in classical case. The inset of Fig. 2 shows di- 
rectly cos n = (Si ■ Si+s n ) I S 2 , giving the details of the 
QSE enhancement of MSRO between several neighboring 
spins. 

The spin correlation length £ is often used to describe 
the strength of MSRO. Despite the magnitude of cos#„, 
£ always tends to infinity when temperature approaches 
Tc from above, so near Tc £ may not be a parame- 
ter that properly reflects MSRO. Above Tc the evalu- 
ation of £ in our formalism is straightforward from the 
long- wavelength behavior of the spin-correlation function 
C q oc l/(q 2 +£~ 2 ). We found that at fixed T/T c , £ always 
increases as S becomes smaller, in contrast to the SM 
where £ is independent of S. At T — 1.17c, 6 = 4.1, 3.7 
and 2.8 for 5 = 1/2, 1 and 00, again demonstrating the 
QSE enhancement of MSRO from another prospective. 

Now let us analyze the SW excitations. In the stan- 
dard magnetism theories such as the random phase 
approximation]^ and its various modified versions |22j. 
SW exist due to the magnetic long-range order, so its 
spectrum is renormalized to zero at Tc- In our formalism 
SW comes from the short-range spin correlations and the 
long-range order is no longer a prerequisite for its exis- 
tence, so SW spectrum can be finite at Tc- In Fig. 2 we 
plot the calculated SW spectrum obtained from Eq. (J3J| at 
Tc- To demonstrate the S'-dependence of the SW renor- 
malisation, the SW spectrum at T = (the FM case) is 
also plotted with all oj q scaled by S. 

Let us estimate the renormalisation factor in the BCC 
Fe2] where SW modes have been observed above the 
middle of the Brillouin zone along the (110) direction, 
Q = (f fO) (lattice constant a = 1). The SW renor- 
malisation factors loq(Tc)/ <jJq (oJq is ojq at T = 0) 
are 0.86,0.76 and 0.60 for S = 1/2,1 and 00, corre- 
spondingly. Experimentally in the BCC Fe loq(Tc)/ 
ujq(0.3T c ) w 0.84 [2] and the difference between ujq and 
luq(0.3T c ) is about 15% 0], so the overall SW renormal- 
isation factor becomes 0.71, and our result for 5 = 1 is 
close to that. Fig. 2 also indicates that w q for smaller 
spins is less affected at elevated temperatures, implying 
that QSE favors the persistence of SW modes. 

Let us now estimate the influence of dynamic effects 
and obtain the relaxation function F(q, uj). Among vari- 



2 



ous analytical approximations for F(q,ui) the three-pole 
approximation |23| seems to be one of the best and it 
has been successfully applied to the typical Heisenberg 
system with large spin S = 7/2 [24L l25j. In this approx- 
imation F(q,ui) is expressed in terms of 81 = (ui 2 ) q and 
<5| = (cj 4 )q/(5j r — 81, where (ui n ) q are frequency moments 
of F depending on the static correlation. The evaluation 
of (w 2 )q is straightforward 23]. (uj 4 ) q oc {[S*,iSt q ]) 
|26l | contains four-spin correlation terms which have to 
be properly decoupled as a product of two-spin corre- 
lations. In the literature the conventional decoupling 
(^i SfSi^Sj } — > Ci-jCf^ and {Sf S m S 3 ) — > Ci m Ci + 
CijCim, ap pro priate for large S, have been applied to ob- 
tain (w 4 )q |23l l2^| . For small S, the spin kernel effect, 
which is neglected in this decoupling, becomes important. 
This QSE can be clearly seen in S = 1/2 case, where for 
i = I or m = j the left side of the decoupled equa- 
tion vanishes while the right side is finite. To take into 
account this QSE we introduce the following decoupling 
procedure 

({Si~, Sf }{S* t , Sj }} — > fafmfiij^lm f° r Ril < Ri m , Rlj , 

I C+ C+ C* — \ . rs rs [ rs rs fi ifS PS r~i 

\°i J l °rn J j I J UJ mjYJ ijJ Im^im^lj ~r J irnJ Ij^ij^lm 

+{$t 3 $i m + 5 im 6i :j )C* l z }, (6) 

where = 1 — 8n/2S. If i,l,m and j are four differ- 
ent sites then Eq.© is the same as in the conventional 
decoupling. QSE occurs when two or more out of these 
four sites are the same. In this case Eq.ljHJ at S = 1/2 is 
exact and is reduced to the conventional decoupling for 
S — » oo. With these results for two opposite limits of 
S and the introduced earlier quantum correction in ff t 
~ l/S, one can expect that Eq.10 will be a reasonable 
interpolation for arbitrary S. By applying this decou- 
pling procedure one can obtain (w 4 )q = (w 4 )q°' + {w i ) q 1 \ 
where (w 4 )^ corresponds to the conventional decoupling 
[23ll2^ | while (w 4 }^ is the quantum correction given by 

^ = 4S^ { iv" ^ [jk{Agl ~ 69k9q+k + 2g «+ k) 
-C k (h k - h k+q )(13J k - 7J q+h )} - (Ug - 9g q )(h - h q ) 
+S- 1 z n JlC n {\ - 7 *)(5C„ + 7C - 65)}, (7) 

n 

where \q is the q-dependent susceptibility, g k — 

En z nJ n C n Jn, and h k = En ^n^nln- 

As a function of u> the relaxation function F(q, oj) has 
either one maximum at uo — 0, if <5f > 281, or three 
maxima at ui — and ui = ±w™ ax , if 8\ < 25f. The latter 



case is often referred as the SW peak at w™ ax 0, HH • 
With such a definition the criteria of the SW existence 
for given q is 8\j8l < 2. Usually w™ ax is slightly larger 
than w q . In the literature the SW peak was also defined 
as (oj) q H3] which is slightly smaller than uj q . Near the 
critical value 8\j8l < 2, the maximum of F at u;™ ax is 
broad. When <5|/<5^ is decreased, the SW peak is more 
pronounced. In Fig. 3 we plot the magnitude of S^/Sf 
for different S as a function of q at Tc- At fixed q, 
6%/6i is always decreased if S becomes smaller. Along 
the (qqO) direction, the critical values of q, when 8\j8l = 
2 , are q CT w 0.30tt, 0.51tt and 0.61tt for S = 1/2,1 and 
oo, correspondingly. Our value of q cr for S = 1 agrees 
with the experiment result in BCC Fe, where SW modes 
above Tc exist only above q ~ ir/2 in (110) direction 
(Fig. 2ofRef. 0). The SW peaks were also obtained in 
the SM0 (with spin independent 8% /Si), but the value of 
<7cr there is considerably higher. Our calculations indicate 
that this theory, which correspond to S = oo, will be 
applicable if QSE is properly taken into account. At 
q = (f §0) for S = 1/2, 1 and oo the ratio <5f/<5? is 
approximately 0.93, 2.2 and 3.6, which are respectively 
well below, close to, and well above the critical value 
82 / 81 = 2. The corresponding dynamic structure factor 
S(q,uj) = lu(1 — e~ UJ / T )~ 1 x q F(q,uj) as a function of uj 
is shown in the inset of Fig. 3. It is clear that at this q 
the well-defined SW exist in the case of S = 1/2, the 
tendency of SW appears for S = 1, and there is no SW 
signal at all for S = 00. Fig. 3 shows that QSE favors 
the persistence of SW with increasing impact for smaller 
spins. In many real magnets S is not large ( S ~ 1 in 
BCC Fe and S w 0.3 in FCC Ni) and we believe that QSE 
plays an important role in the MSRO and the magnetic 
excitations above Tc, especially in the itinerant magnets. 

In conclusion, we analytically demonstrated the pres- 
ence of MSRO in the Heisenberg model and identified 
the importance of quantum spin effect on MSRO for 
ferromagnets above Tc- By extending the second-order 
Green's function technique to arbitrary S we found that 
for a system of small spins the quantum effects greatly 
contribute to the MSRO and enhance its influence. The 
spin dynamics investigation developed from the conven- 
tional method of moments further confirms that QSE fa- 
vors the persistence of spin wave excitations. We demon- 
strated that this previously neglected QSE removes the 
long-standing controversy between theory and experi- 
ment regarding the presence of MSRO and SW in Fe 
and Ni above Tq and clearly indicates that the current 
prevailing point of view of finite temperature magnetism 
should be reconsidered to properly include MSRO and 
quantum effects. 
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Fig.l. T C /T^ F and E C /E Q as a function of S from 
SM (lines), the HTE methods (close symbols) and our 
formalism (open symbols) in BCC structure in the NN 
coupling case. 

Fig. 2. The calculated SW spectrum u q for the differ- 
ent S at T c . The dashed line is w q at T = (FM case). 
The inset shows cos 9 n from nearest to fifth-nearest neigh- 
bors. 

Fig. 3. The calculated S%/5f for the different S at Tc as 
a function of q. The criteria of SW <5f /<5i = 2 is marked 
by the dashed line. The inset shows S(q, u>) at q = (5 -f 0). 
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